Critical behavior of charmonia across the phase transition: A QCD sum rule approach 



Kenji Morita^B and Su Houng Lee 1 '!!] 

institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 

(Dated: May 14, 2008) 

We investigate medium-induced change of mass and width of J/tp and r\ c across the phase transi- 
tion in hot gluonic matter using QCD sum rules. In the QCD sum rule approach, the medium effect 
on heavy quarkonia is induced by the change of both scalar and twist-2 gluon condensates, whose 
temperature dependences are extracted from the lattice calculations of energy density and pressure. 
Although the stability of the operator product expansion side seems to break down at T > 1.06T C 
for the vector channel and T > 1.04T C for the pseudoscalar channel, we find a sudden change of the 
spectral property across the critical temperature T c , which originates from an equally rapid change 
of the scalar gluon condensate characterized by e — 3p. By parameterizing the ground state of the 
spectral density by the Breit-Wigner form, we find that for both J/tp and rj c , the masses suddenly 
decrease maximally by a few hundreds of MeV and the widths broaden to ~ 100 MeV slightly 
above T c . Implications for recent and future heavy ion experiments are discussed. We also carry 
out a similar analysis for charmonia in nuclear matter, which could serve as a testing ground for 
observing the precursor phenomena of the QCD phase transition. We finally discuss the possibility 
of observing the mass shift at nuclear matter at the FAIR project at GSI. 

PACS numbers: 14.40.Gx,11.55.Hx,12.38.Mh,24.85.+p 



I. INTRODUCTION 



In-medium change of spectral properties of heavy 
quarkonia is one of the interesting problems in recent 
hadron physics. Firstly, the recent relativistic heavy ion 
collision experiment at the Relativistic Heavy Ion Col- 
lider (RHIC) reveals exciting nature of the QCD matter 
through a number of observations [E 0, H, EL III • However, 
there are many open questions in both experimental facts 
and theoretical understandings of QCD matter. Hence, 
it is important to establish appropriate experimental ob- 
servables that reflect consequences of deeper theoretical 
understanding of the matter. Heavy quarkonia have been 
regarded as one of the most suitable diagnostic tools in 
this respect, since the suppression of J/tp yields would 
reflect the Debye screening phenomenon caused by the 
deconfinement phenomenon in the quark-gluon plasma 
(QGP), as was originally argued by Matsui and Satz 
Q. Until now, quarkonium production, especially that 
of J ftp, in relativistic heavy ion collisions have been ex- 
tensively studied both experimentally @, and theoreti- 
cally [3, US- However, a remarkable progress comes from 
recent lattice QCD calculations, which indicate that con- 
trary to the earlier expectation the J/ip will survive as 
a bound state even in the QGP up to T ~ 1.6 — 2T C 
[Til H3 , [l3| , which was anticipated before based on the 
non-perturbative nature of QGP [13] • Nowadays, the 
state of matter at this temperature region has been char- 
acterized as "strongly coupled" QGP (sQGP). Hence, 
there will be change of spectral properties even for heavy 
quark system which has to be considered in interpreting 
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experimental observables. 

Secondly, charmonium in a nuclear medium is also 
an interesting issue. In relativistic heavy ion collisions, 
we need knowledge of quarkonium-nucleon interaction to 
discriminate the suppression by QGP from the "cold nu- 
clear matter effect" induced by such an interaction. Fur- 
thermore, multi-gluon exchange can lead to an attractive 
interaction between cc and a nucleon, which may result 
in a bound state of charmonium and a light nuclei, as 
pointed out by Brodsky et al. 15]. It should be also noted 
that the Panda experiment at GSI-FAIR plans reaction 
of anti- protons with nucleus target, which will yield char- 
monia in the nuclear matter. It could serve as a testing 
ground for observing the precursor phenomenon of the 
QCD phase transition. 

In this paper, we investigate change of mass and width 
of J/tp and rj c induced by strongly interacting hot glu- 
onic matter and by nuclear medium using QCD sum rule. 
QCD sum rule provides a systematic procedure for study- 
ing hadrons from a viewpoint of the asymptotic freedom 
in QCD [l6l.[l7|. Since the QCD sum rule can take non- 
perturbative effects into account through the condensate 
terms, it is a suitable theoretical tool of the current study. 
Indeed, QGP at T < 3T C cannot be understood using 
perturbation theory alone [lq ]. Furthermore, the sum 
rule is more promising for heavy quark systems because 
we do not have to take the quark-antiquark condensate 
into account unlike light quark systems. In this respect, 
the sum rule has been applied to charmonium and bot- 
tonium. Shifman et al. established the framework in 
Ref. [lfj, [l7| and Reinders et al. extended it to deep Eu- 
clidean region Q 2 = —q 2 > [l9(, in the case of vacuum. 
As for the quarkonia in-medium, One of us together with 
Furnstahl and Hatsuda have investigated the mass shift 
of J/tp in hot hadronic matter [2(J, using a QCD sum 
rule approach, where the temperature effect was intro- 



2 



duced to the perturbative Wilson coefficient through the 
scattering terms. A consistent formalism at lower den- 
sity was developed by one of us [2l| and independently 
by Hayashigaki [12] to study the mass shift of Jji\> in 
nuclear matter. 

Along this direction, we investigated the mass shift and 
width broadening of J ftp in hot gluonic plasma (GP) [23| 
just above the phase transition by consistently using the 
exact temperature dependencies of condensates from lat- 
tice calculation. In the present paper, as a subsequent 
paper to Ref. (23[, we present details of the analysis, fur- 
ther application to r\ c and to spectral changes in nuclear 
matter. 

In the next section, we will give an explanation of the 
QCD sum rule for heavy quarkonium in medium used 
in the present work. Section IIIII and IIVI describe the 
details of the numerical computations of the sum rule 
for hot gluonic matter and nuclear medium, respectively. 
Section [V] is devoted to discussion and summary. 



II. QCD SUM RULE FOR HEAVY 
QUARKONIUM 

In this section, first we review the sum rule for heavy 
quarkonium in vacuum [191 ]. Then we introduce the ex- 
tension to finite temperature and nuclear medium cases, 
in which medium effect is eventually induced only by the 
expectation values of gluonic operators without any addi- 
tional change in the operator product expansion (OPE). 



A. Moment sum rule in vacuum 

We start with the time-ordered current-current corre- 
lation function for J channel 



n J (g)^ / d 4 xe^<T[/(*)/(0)]), 



(1) 



where we consider J = P (pseudoscalar) and V (vector) 
current of the heavy quark. Namely, j p — icj^c and 
jY — c-f^c, for charm. The expectation value (• • ■ ) is 
taken for the vacuum. If we go to deep Euclidean region 
Q 2 = —q 2 » 0, the product of the current can be ex- 
panded via operator production expansion (OPE) [24| . If 
we denote U(q 2 ) such that W(q) = (qf q v - q 2 g^ v )i\{q 2 ) 
for the vector current, n(g 2 ) can be written as 



irV) = E c «<°«> 



(2) 



where O n are the operators of mass dimension n renor- 
malized at scale /i 2 and C„ are the Wilson coefficient. By 
virtue of much heavier quark mass than the confinement 
scale, heavy quark operators, such as m c cc for dimension 
4, are rewritten in terms of gluonic operator with a factor 
of l/m c via heavy quark expansion [l6t [25j. Hence, only 



gluonic operators contribute to the OPE for the heavy 
quark currents. 

On the other hand, the correlation function ^ is re- 
lated to its imaginary part through the dispersion rela- 
tion 



Imn J (s) 



ds 



(3) 



where we ignore +is in the denominator of the inte- 
grand since q 2 = —Q 2 < 0. Taking n times derivative 
of Eqs. (J21) and ([3]) as 



M„ J (Q 2 ) 



it 



d 

dq^ 



we obtain the n-th order moment for the OPE side 

M„ j (Q 2 )ope - 4l(0[l + ai(0a s + b J n ^)cj> h ], 
and that for the phenomenological (dispersion) side 

1 r°° ImfT 7 (s) 



M„ J (Q 2 )p ho 



(s + Q 2 Y 



+ i 



ds. 



(4) 



(5) 



(6) 



Here, we have introduced a dimensionless scale variable 
t = Q 2 /im 2 c . In Eq. ©, A J n (0, and b J n (0 are 

the Wilson coefficients which correspond to bare loop di- 
agrams, perturbative radiative correction up to order a s , 
and scalar gluon condensate, respectively. These coeffi- 
cients were derived in Ref. [191 ] and we summarize them 
in the Appendix. 

In evaluation of spectral properties, we take the ratio 
of the (n— l)-th moment to the n-th moment and equate 
the OPE side with the phenomenological side. Then we 
obtain the sum rule 



Mi 



OPE 



Mi 



(7) 



phcn. 



which relates the hadron properties (r.h.s.) with asymp- 
totically free QCD (l.h.s.) 



B. Moment sum rule for the hot gluonic medium 

In this paper, we firstly consider the gluonic medium 
at finite temperature around T c . Then, the expectation 
value in Eq. Q is taken as (O) = Tr( e - /5H C')/Tr(e~ /3ff ). 
Hereafter, we set both medium and cc at rest. We de- 
note q^ = (uj, q) and take q — > limit. In this case, the 
transverse and the longitudinal components of the corre- 
lation function for the vector channel are simply related 
with n T = u 2 n L and n L = n^/(-3w 2 ). We denote the 
longitudinal component as Ii J (w) for the vector channel. 

At finite temperature, retarded correlation function is 
related to the spectral function 26]. In the Euclidean 
region uj 2 < 0, the retarded correlation function U r (uj) 
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becomes n(w 2 ) and the dispersion relation is given by 

[HI13 

n-v)= r d u 2 -^, (8) 

J - U Z - LO l 

where p(u) is the spectral function connecting with the 
imaginary part as 

p( U Hitanh(^)lmITV). (9) 

Then Eq. ((8]) reduces to the vacuum case [Eq. ((3|)] when 
Imfi J (it 2 ) has nonzero value only at u > T. Since we 
are interested in charmonia for which the mass is much 
larger than temperature considered here, this condition 
seems to be appropriate one. However, there are formally 
two additional terms in the finite temperature spectral 
function (28j . One is the continuum part which also ex- 
ists in the case of vacuum. Following the prescription in 
Ref. , we can suppress contribution from this part as 
described later because this part has finite values beyond 
some threshold. The other part arising from scattering 
of the current with quarks in medium is proportional 
to S(u 2 ) and the contribution grows up with T in the 
hadronic medium [20l ]. However, since we are considering 
the gluonic medium in which there are no (anti-)quarks 
which annihilate with the current, such a scattering term 
does not appear. Hence, we can use the same expres- 
sion of the phenomenological side with the vacuum case 
[Eq. ([5])] for charmonia in the hot gluonic medium. 

As for the OPE side, there is an important change 
from the vacuum to the medium case. Since we have no 
longer Lorentz invariance, non-scalar operators have non- 
vanishing value [27} ■ In the present case, twist-2 gluon 
operator has leading contribution and the n-th order mo- 
ment of the OPE side [Eq. |(SJ)] should be modified to 

m„ j (q 2 )ope = Aim + an(£K + rtiOfa + ci(OU 

(10) 

where c n and 4> c are the Wilson coefficients and the 
medium expectation value for the twist-2 operator. Since 
we are considering the heavy quark systems, only the 
condensate terms are temperature dependent as long as 
T <C m c , \Q\ [13, Hence, the Wilson coefficients are 
the same as in the vacuum case. In the following, we show 
that the gluon condensates 0b, c are written in terms of 
thermodynamic quantities which can be extracted from 
lattice QCD data. 

If we define these condensate terms as 

Gq(T) = l—G a G a,J,u \ , (11) 

\ 7T K IT 

(n^u" - \g^ v J G 2 (T) = (?!lGWG wp } T , (12) 
where is the 4-velocity of the medium and taken to 
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FIG. 1: (Color online) Temperature dependent coupling con- 
stants extracted from lattice QCD. The boxes denote the lat- 
tice data points of a qq (r scrccn , T) taken from Ref. [29]. The 
solid line is drawn by Bezier interpolation of the lattice data 
points. The dotted line shows the case of Eq. (|20p . 

be m m = (1,0,0,0), explicit forms of </>b, c are given as 

4tt 2 

0b = 9(4^F Go(T) ' (13) 

^ = ^7l4^ G 2(T). (14) 
3(4m^) z 

Actually it is possible to calculate the condensates (fTTj) 
and (fT2")) directly using lattice QCD, but we do not 
adopt such an approach here. The gluon condensates 
generally consist of the perturbative piece and the non- 
perturbative piece. At zero temperature, the condensate 
term appearing in QCD sum rules is the non-perturbative 
piece only and it is shown that the non-perturbative part 
extracted from lattice QCD by subtracting the pertur- 
bative part is indeed consistent with the value of the 
condensate determined from QCD sum rules for charmo- 
nium [13, HH, H3] • Similar consideration holds also for the 
finite temperature case [33], in which we would have to 
subtract out the perturbative part at T ^ if we directly 
calculated the non-perturbative condensates from lattice 
QCD. In this paper, since we are putting all the tem- 
perature dependencies in the condensates, including the 
perturbative and non-perturbative contributions, we can 
just extract total temperature dependencies of the opera- 
tors from the lattice. This is possible by noting that the 
scalar gluon condensate and twist-2 gluon condensates 
are respectively just the trace part and symmetric trace- 
less part of the energy momentum tensor. This energy 
momentum tensor is well calculated on the lattice from 
the pressure and energy density of the plasma through 
the following equation, 

T aP = (e+p) (u a u? - \g a ^j + \{e - 3p)g^. (15) 

The scalar condensate can be related to the trace part 
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through the trace anomaly as 



use the 2-loop perturbative running coupling form 



rru_/m 



(16) 



with /3(g) being the beta function, (3(g) = 

— 4§^3- (33 — 2Nf) for 1-loop, N{ flavors, and N c colors. 
Using the above expression with N{ = and N c — 3 for 
the beta function and recalling that = e — 3p, we 
obtain 



G Q (T)=Gr--(e-3p) 



(17) 



where Gg ac is the value of the scalar gluon condensate 
in vacuum [34j. As for the twist- 2 part, the symmet- 
ric traceless part of the energy-momentum tensor is the 
gluon operator 



(18) 



Hence we can identify the traceless part of the energy 
momentum tensor to (e + p) as given in Eq. (|15p . ^From 
Eq. (TT2"!) . the twist-2 part becomes 



G 2 (T) 



a s (T) 



(e + p), 



(19) 



so that G 2 (T) is proportional to the entropy density of 
the system s = (e +p)/T. We extract the temperature 
dependent quantities e, p (3f| and a s (T) [2§| from lattice 
calculations for the pure SU(3) system. In order to con- 
struct G2, we need the temperature dependent effective 
coupling constant. The coupling constant, however, can- 
not be uniquely determined by lattice QCD [2^| . Ref. 
presented four kinds of the coupling constant extracted 
from the color singlet heavy quark-antiquark free energy. 
Two of them are measured in the short distant regime 
and the others are done in the long distant regime. In the 
former, one is from the free energy and the other is from 
the spatial derivative of the free energy (force). Both cou- 
pling constants are almost independent of temperature at 
short distance, r < 0.1 fm. While the former goes to neg- 
ative value at larger distance due to the remnant of the 
confinement force, the latter shows temperature depen- 
dent maximum value, at which the distance is denoted 
by hereon- Here, we adopt the latter one, a qq (r, T) at 
f = ^screen as one of reasonable coupling constants since 
it characterizes the relevant length scale for the separa- 
tion of short distance regime from long distance one. On 
the other hand, the long distant regime is based on a 
fit of the free energy to the Debye-screened functional 
form which has two coupling parameters, Coulomb force 
strength a(T) and screening a(T). Although both of the 
coupling constants show reasonable temperature depen- 
dencies and agree each other at T > 6T C , we adopt a(T) 
because the Coulomb force strength is not relevant for 
characterizing the long distance non-perturbative physics 
at temperature considered here. Unlike a qq , the uncer- 
tainty in the result of a(T) is too large. Therefore, we 
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(20) 

j^jg _ 1.14 and rescale this as &[T) = 
2.095a po ' rt (T) [11]. Here we put T c = 264 MeV The 
two coupling constants as a function of temperature are 
displayed in Fig. [T] As explained later, we will consider 
only temperature region near T c in this paper. Hence, 
a qq is stronger than a(T) throughout analyses in this 
paper. 




FIG. 2: (Color online) Gluon condensates near T c . 

The resultant gluon condensates Go and G 2 for two 
cases of the coupling constant are shown in Fig. [Sf . For 
G , we use G^ ac = (0.35GeV) 4 ~ 0.015GeV 4 . We can see 
that Go decreases as temperature increases and reaches 
less than half of the vacuum value at T/T c ~ 1.04. It be- 
comes negative at higher temperature but remains posi- 
tive in the temperature region considered here [331 ] . 



C. Moment sum rule for the nuclear medium 

In this case, the medium consist of nucleons, thus we 
do not have to worry about the scattering term. As far 
as we follow the same method to suppress the contribu- 
tion from the continuum, we can use the same form of 
the phenomenological side with the vacuum and finite 
temperature cases. 

Thus, since the medium effect is similarly imposed on 
the gluon condensates, difference in the nuclear matter 
case from the case of hot gluonic matter is in the explicit 
form of 0b c- In order to evaluate the expectation value 



1 We have renewed the extraction from lattice data by improving 
the resolution, so that the present values are slightly different 
from those of Ref. [H. 
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FIG. 3: (Color online) Gluon condensates in nuclear matter. 
Thick solid and dashed line show the scalar and twist-2 con- 
densates as a function of density normalized by the normal 
nuclear density. Thin lines are the finite temperature case at 
T = T c for a comparison. 



The form of 4>b is the same as the hot gluonic matter case 
but now the expectation value is taken through Eq. (|22p . 
We depict the density dependence of the gluon conden- 
sates based on Eqs. and (gSJ) in Fig. [3] The twist-2 
case is re-normalized so that it corresponds to the finite 
temperature case (fl"3)l . We can see that the change of 
the scalar condensate reaches as large as T — T c case at 
p ~ 5po but is much smaller at the normal nuclear den- 
sity. The twist-2 contribution is much smaller than that 
of the finite temperature case. 



D. Phenomenological side 

In the phenomenological side, we use a simple pre- 
scription which describes the lowest lying resonance 
in each channel. For charmonium, previous studies 
pl M M HE m focused on the mass and ignored the 
small but finite width of J/ip and rj c . In this case, the 
imaginary part of the polarization function in Eq. ([B]) is 
simply parametrized by 



for the ground state of the nuclear matter, we employ the 
linear density approximation (36J: 



(O)n.m. = (O} + 



Pn 
2m at 



(N\0\N), 



(21) 



where pn and mjv are the normal nuclear matter density 
and the nucleon mass, respectively. The nucleon state 
\N) is normalized as (N(p')\N(p)) = 2p {27r) 3 S 3 (p - p'). 
Then, the scalar condensate becomes 21il 



7T / n.m. 



7T p / 



-m° N p N (22) 



where m a N ~ 750 MeV is the nucleon mass in the chiral 



limit [37[. The traceless and symmetric twist-2 operator 
is given as plf . 



N(p) 



< - r a<7 ( - r 

7T 



after 



N(p)) = -{p a p /} -\ 9 yj ^a g 

""(23) 

where Aq is related to the moment of the gluon distri- 
bution function 



A G (p 2 ) = 2 / 
Jo 



dx xG(x, p 2 ). 



(24) 



Following Ref. [H|, we take A G (8m 2 c ) ~ 0.9. While G 2 
at finite temperature is related to the entropy, this corre- 
spondence does not hold in the nuclear matter case. Note 
that Eq. (p~8|) does not contain the quark sector. Using 
these expressions, the condensate terms which appear in 
Eq. © finally result in [H| 



in 2 
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7T 'n.m. 
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2 m NPN- 



(25) 
(26) 



ImII(s) = Jq5(s — to 2 ) + corrections, 



(27) 



where we ignore the channel subscript J. This spectral 
function immediately leads to the moment 



fo 



ir(m 2 + Q 2 ) 



2\n+ 



T [l + S n (0}. (28) 



The correction term in Eq. [f2"7]) is absorbed in 5 n (£). By 
taking the ratio as in Eq. ([7]) , we can remove the constant 
/o from the equation. To obtain the mass of lowest lying 
resonance, we need to choose sufficiently large n such 
that (1 + <J n _i(£))/(l + <5„(0) is close to unity. Then 
the ratio does not depend on the details of the correction 
term which contains higher resonances and continuum, 
and the mass is simply given by 



M n _i(Q 



(29) 



Previous analyses rely on this formula. 

In this work, we extend the above formulation to in- 
clude finite width. Here, we employ the simple relativistic 
Breit-Wigner form 



ImII(s) 



(s - to 2 ) 2 + s r 2 



corrections. 



(30) 



As in the r = case, we can eliminate the unnecessary 
constant and the effects of the correction term by taking 
the ratio of the moment and choosing appropriately large 
n. In the practical analyses of the sum rule, our task is to 
find values of (m, V) which satisfy the sum rule [Eq. ([7}]. 
Generally there are infinite numbers of the pairs of (to, T) 
because the sum rule provides one equation with respect 
to the two quantities which we want to know. Hence, 
without additional constraints, the sum rule can provide 
only relation between to and T as in the case of light 
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vector mesons [38J. Here, before the practical calculation, 
we discuss the relation between the mass and the width 
which comes from the phenomenological side, Eq. ([30]) . 

In calculation of the moment ratio of the phenomeno- 
logical side, we need to compute the dispersion integral in 
Eq. ([S]) with the spectral function in Eq. ([50]) . Since the 
width of the ground state charmonium is much smaller 
than its mass, we need careful treatment in numerical 
integration. To achieve good accuracy, we performed 
Monte-Carlo integration based on the VEGAS algorithm 
39]. In our calculation, typical relative numerical uncer- 
tainty evaluated from the standard manner in the Monte- 
Carlo integration is order of for 10 6 events with 
m = 3 GcV and T = 1 MeV. As expected, this accuracy 
becomes better as T increases. 
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FIG. 4: (Color online) Moment ratio of the phenomenological 
side as a function of V. Upper panels are for £ = 1. Left and 
right panels denote the case of n = 9 and n = 14, respectively. 
Lower ones are for £ = 3 with n = 14 (left) and n = 19 (right). 



We plot the T dependence of the moment ratio for var- 
ious mass values from 2.9 GeV to 3.1GeV in Fig. 2] Here 
we show the result for two values of £, £ = 1 and 3. In 
each £ value, we choose two values of n, of which are the 
typical values for the analyses below, to see n dependence 
of the moment ratio. First, comparing the left (smaller 
n) side to the right (large n) side, we can see that T 
dependence of the moment ratio becomes stronger as n 
increases. As we will see later, larger n is suitable for 
evaluating mass at higher temperature. Hence, this fact 
means that, as the temperature increase, the system be- 
comes more sensitive to the width. Second, the moment 
ratio decreases monotonically as the width increases if 
mass is unchanged. It also decreases as the mass de- 
creases but the width dependence is much weaker. For 
instance, let us suppose that we obtain 1 GeV 2 decrease 
of the moment ratio from the OPE side for £ = 1. If 
mass stays constant, the width must broaden to larger 



than 100 MeV while it corresponds to about 100 MeV 
mass reduction in the case that the width remains in its 
vacuum value. Finally, as is shown in comparison of the 
upper-right with the lower- left, the width dependence be- 
comes weaker if we choose larger £. Its consequence will 
be discussed in the next section. 



III. CHARMONIUM IN HOT GLUONIC 
MATTER 

In this section, we present the result of the analysis for 
the hot gluonic matter. The parameters of the theory are 
a s and m c . Hereafter, they are set to 0.21 and 1.24 GeV 
at £ = 1, that are taken from [2l| , respectively. 

We begin with fixing n such that the moment ra- 
tio of the OPE side takes its minimum value for each 
temperature. As briefly mentioned before, we need to 
choose moderately large n so that contribution from ex- 
cited states and continuum can be sufficiently suppressed. 
Therefore, this ratio should approach a constant value at 
adequately large n. However, in the OPE side contribu- 
tion from higher dimensional operators will be important 
at large n. As such n value that the moment ratio takes 
its minimum value, pole dominance and truncation of the 
OPE are valid and the ratio is close to the real asymptotic 
value, as have been extensively studied in the vacuum 
case [19j . 

We display the moment ratio for the OPE side 
[Eq. (Unj)] in Figs. O and [5] with the gluon condensates 
shown in Fig. [5] 

Figure [5] shows the moment ratio for the vector chan- 
nel. The left and right column show the case in which we 
use a qq and a, respectively. Comparing different £ cases, 
we can see that the stability of the moment becomes bet- 
ter as £ increases. But the values of n which give the 
stability to the moment ratio also becomes larger. As 
previously reported in [23[ , the stability is only achieved 
near T c and the stronger coupling, which is a qq in this 
temperature region, gives worse stability. By increas- 
ing £, we can improve the stability a little. While it is 
achieved only up to 1.04T C for £ = 0, the moment ratio 
remains stable up to 1.06T C for £ = 3. 

We can see the similar situation in the pseudoscalar 
channel depicted in Fig. [S] However, the moment ratio 
is less stable than the vector case. In the pseudoscalar 
case, even the best case (using a and £ = 3) can stabilize 
the moment ratio only up to 1.04T C . 

Note that the lack of stability does not necessarily 
mean dissociation of the charmonia. The reason of such 
instability can be clearly seen in the each terms of the 
OPE [Eq. (flU)) ]. of which each term must be much less 
than unity for convergence. These terms are displayed 
in Figs. [Tim We can see that all the coefficients grow up 
with n. An important feature in all the coefficients is that 
increasing £ clearly keeps their value smaller. Among 
these three, only c„(£)0 c always has positive sign and 
its magnitude increases with temperature. These two 
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FIG. 5: (Color online) Moment ratio for the OPE side for 
the vector channel (J/ip). Each panels show different £ and 
coupling constant case. The symbols stand for different tem- 
perature. 



FIG. 6: (Color online) Same as Fig. [5] but for the pseudoscalar 
channel (rj c ). 



features are opposite to b n (£)<j)b, in which the sign is al- 
ways negative and the value seems to approach to zero as 
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FIG. 7: (Color online) Radiative correction term a„(£)a s in 
the OPE. Upper panel shows the vector case and lower one 
shows the pseudoscalar case. 



temperature increases. In comparing the two channels, 
one finds that there are no significant differences. Hence, 
the stability will be determined by a delicate balance be- 
tween coefficients and its breakdown will be caused by 
rapid increase of c„(£)0 c . 

Now we proceed to the determination of mass and 
width. The values of n are listed in Table |TJ Note that 
the stability achieved at the highest temperature is am- 
biguous in some cases; for example, J/ip of £ = 2 with 
a m case is stable at T/T c = 1.05 with n = 22. However, 
as seen in Fig. [5l the moment ratio is almost constant 
in such large n region and never rises up as lower tem- 
perature cases do. Such a vague stability is also seen in 
other cases. Hence, we note that mass and width values 
evaluated on the basis of such a stability are less reliable 
in the analyses below. 

Once n and £ are fixed, we can compute the mass and 
the width by making use of Eq. ([7]) . For a fixed moment 
ratio of the OPE side, we firstly compute the mass in the 
limit of r — > using Eq. (|2"5|) . By virtue of the monotonic 
behavior of the moment ratio of the phenomenological 
side shown in Fig. [4] we can safely calculate the mass in 
the case of finite width by numerically solving Eq. (JT]) 
with Eq. ([501. 

We plot the relation between the mass shift and the 
width at various temperatures in Figs. [TU] and [TTJ We 
can see the almost linear behavior of the width as a func- 
tion of the mass shift. Note that the vacuum mass dif- 
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FIG. 8: (Color online) Scalar condensate term &n(£)</>b- Left 
and right column stand for the vector and the pseudoscalar 
case, respectively. Symbols denote different temperature 
cases. 



fers for different £. We do not perform fine tuning of 
the parameters so that the real vacuum mass is repro- 
duced. Although there are some exceptions for the linear 
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relation, especially small £ and high temperature cases, 
these come from the vague stability we mentioned be- 
fore. Hence, we can conclude that the mass shift and the 
width have the linear relationship as far as QCD sum 
rules properly work. The other important aspect is tern- 
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perature dependence of the mass shift and the width. 
We cannot know how the mass and the width behave in 
the real situation, since we cannot simultaneously deter- 
mine both of the mass and the width within the current 
framework only. Here, we investigate two extreme cases; 
r — > limit and 8m — > limit. 

The results are shown in Figs. [12] and HH In these 
figures, we plot the results of T > 0.9T C . Figure [T2l 
shows the remarkable behavior of the mass shift; The 
mass does not change up to T ~ T c but it suddenly be- 
gins to decrease across T c . This fact clearly reflects the 
temperature dependence of the gluon condensates which 
represent the phase transition. Above T Cl the mass de- 
creases with temperature almost linearly. This feature 
is common for both J/t/j and r\ c . Though small £ re- 
sults, especially £ = 0, show more rapid decrease, the 
curves become almost parallel among large £ results, as 
a consequence of the better stability. From the nature 
of the phenomenological side shown in Fig. 2J this case 
corresponds to the maximum mass shift. The mass shift 
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FIG. 10: (Color online) Relation between mass shift Sm = 
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figure shows different £ case for two cases of the coupling 
constant, a qq and a. 



shows ~ 50 MeV reduction from vacuum to T c , and it 
increases additionally by ~ 20-50 MeV as temperature 
rises by 0.01T C . Consequently , it becomes 100-300 MeV 
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at T = 1.04T C . 

Similarly, Fig. Q]|] shows that the width begins to in- 
crease with temperature across T c if no mass shift takes 
place. This also shows almost linear dependence on tem- 
perature above T c . Though some exceptions can be seen 
in the small £ results, which are also indicated in Figs.fTUl 
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Mass shift of JA|/ and r| c , for r=0 
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FIG. 12: (Color online) Temperature dependence of the 
masses in the F — > limit. Symbols stand for different £ 
values. 



and ITT} these behaviors come from the vague stability 
appearing as too large n in Table UJ Hence, we can con- 
clude that the width increases linearly with temperature 
above T c if the mass remains unchanged. Since we did 
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not do fine tuning of the parameters for each £, the values 
of mass shift and width differ for different £. However, 
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the qualitative features do not depend on £ where the 
stability is reliable. This shows robustness of our analy- 
sis. A realistic change at each temperature should be a 
combined decrease in mass and increase in width, whose 
values are smaller than their maximal changes obtained 
here. However, to determine the realistic combination, 
we need to have an additional constraint between the 
changes in the width and the mass, or input the thermal 
width from another calculation. From Fig. [2l one may 
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FIG. 14: (Color online) 5m of J/ip without twist-2 term Gi. 
See text for detail. 

think that the dominant contribution to the change of 
mass and/or width is the scalar gluon condensate which 
exhibits sudden decrease around T c . However, G2 has 
also similar behavior since it relates to the entropy den- 
sity. Though the value of G2 around T c is smaller than 
Go because of prefactor a s /Tt, the relative contribution 
to the moment becomes larger as T increases. In order 
to see the contribution clearly, We show the mass shift 
of J/tp without G2 term for £ = 1 together with the two 
different coupling cases in Fig. [T3J We can see that al- 
most half of the mass shift is caused by decrease of G2. 
Clearly larger G2 value in which a qq is adopted as cou- 
pling constant leads to larger mass shift. At T — 1.04T C , 
a qq is about 0.1 larger than a. This difference makes 
the mass shift 30 MeV larger in the £ = 1 case. Unfor- 
tunately present analysis is limited to the temperature 
region around T c , the role of twist-2 term will become 
more important at higher temperature. 



of Eqs. (fl~3]) and (fl~4| . respectively. With the common 
parameter set, the condensates are <fib — 1-74 x 10~ 3 
for vacuum, 1.64 x 10~ 3 for the nuclear matter, and 
4> c = -1.28 x 10~ 5 . 

As previously shown in Ref. [2l[ , the change of mass, 
which is identical to the change of the moment ratio of the 
OPE side [Eq. (f2T))) ]. is not as large as in the hot gluonic 
matter case. Thus we do not have to worry about the 
stability of the OPE. Nevertheless, increasing £ improves 
the validity of the OPE. We will show the results for < 
£ < 3 as well as in the hot gluonic matter case to show 
the robustness and the consistency of the calculation. 



TABLE II: List of n which stabilize the moment ratio for the 
nuclear matter 
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IV. NUCLEAR MATTER 



FIG. 15: (Color online) Relation between mass shift and 
width in the nuclear matter. 



In this section, we analyze change of mass and width 
of the charmonium induced by nuclear medium with 
the same framework that was implemented in the pre- 
vious section. Here, we use Eqs. (|2l))) and (|26p instead 



See Ref. [40H for a recent investigation. 



We list values of n in which the moment ratio of the 
OPE side for nuclear medium becomes minimum in Ta- 
bic mi The vacuum case has been already shown in Ta- 
bic |TJ Comparing these two cases, we can see that the 
values of n are the same except for a few exception in 
the J/ip case, by virtue of the small shift of the gluon 
condensates in the nuclear matter. In such exceptional 
cases, difference of the values of the moment ratio from 
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the same n value case with the vacuum is almost negli- 
gible, i.e., the moment ratio is almost constant around 
these n. 

We plot the width T as a function of the mass shift 8m 
in Fig. [15] as well as in the GP case. In both J /tp and rj c 
cases, smaller £ than 1.5 show larger mass shift and width 
broadening but larger £ results agree each other. From 
the stability argument, larger £ results will be more reli- 
able. Then, possible mass shifts are maximally -7 MeV 
for J/ip and -4 MeV rj c while maximum widths are 10 
MeV for J/i/j and 6 MeV for rj c . 



V. DISCUSSION AND SUMMARY 



In Sec. IIII1 we have shown that mass decreases sud- 
denly across T c and the shift reaches maximally a few 
hundred MeV above T c in the hot gluonic matter. Al- 
ternatively, width can also maximally broaden to ^200 
MeV. Although our analysis cannot determine both of 
mass and width simultaneously, this is a notable result 
which should be examined in the present and future ex- 
periments. In fact, a next to leading order QCD calcula- 
tion shows that the thermal width of J / i/j slight above T c 
is smaller than a few 10 MeV [ll|,|43|. Hence a large mass 
shift will take place. Note that such a large mass shift has 
been expected from different points view; an AdS/QCD 
analysis shows a sudden drop of mass at the phase tran- 
sition [43[. In Ref. [43j], although the mass begins to 
slowly increase at higher temperature, the temperature 
region investigated in the present paper corresponds to 
the critical region. Sudden reduction of the asymptotic 
value of the potential seen in lattice QCD |44| leads to 
lowering of the bound state energy [45|]. Recent lattice 
QCD calculation based on the maximum entropy method 
also shows survival of the peak in the spectral function 
above T c [46| but the resolution is still insufficient to dis- 
cuss shift and broadening of the peak. Since our results 
access only near T c , we are still far from the complete 
understanding of the behavior of the charmonium in the 
deconfined medium. In the most plausible picture from 
the current understanding, charmonia are melting at very 
high temperature expected in the early stage of the heavy 
ion collisions at RHIC and LHC. Then the pairs of heavy 
quark and antiquark form the bound states at a certain 
temperature which depends on quantum number. The 
temperature is expected as ~ 2T C for J /if) at RHIC [47T |. 
After charmonia are produced, they will dissociate by 
collisions with partons. If this phase lasts long enough 
compared to the inverse of the width, the charmonia can 
decay in the medium. In fact, the lifetime of the partonic 
medium is about 4-5 fm/c in a hydrodynamic calculation 
for the central Au+Au collisions at the maximum RHIC 
energy [48|. This will be much longer at LHC. From 
Fig. [TUl we expect ~ 200 MeV J/i/j mass reduction in 
the case of the small decay width. This shift is larger 
than experimental mass resolutions (~ 35 MeV for di- 
electron channel of PHENIX at RHIC H, 33 MeV for 



dielectron channel and 75 MeV for dimuon channel of 
ALICE at LHC H|). 

Alternatively, statistical hadronization near phase 
boundary has been also examined 50]. In this case, 
the number of produced charmonium will be enhanced 
if the notable mass shift occurs. For example, there 
may be a factor of 2 enhancement for T = 170 MeV 
and 8m = —100 MeV since the enhancement factor is 
given by e~ Sm / T , This enhancement might be observed 
by comparing particle ratio. 

As for the nuclear medium result, we have extended the 
analysis carried out in Ref. [H| to the one which takes 
account of finite width. We have also shown the results 
for different £ values. Since we have given the relation be- 
tween the mass shift and width, we can estimate the mass 
shift in the presence of finite width effect by considering 
the dissociation cross section of the charmonium by nu- 
cleon. Provided the Fermi momentum is p-p ~ 250 MeV 
and the cross section is crj/^_ N ~ 2mb, the decay width 
r = (uj/^-nVi-cIPn) becomes ~ 1.3 MeV for charmo- 
nium at rest. The cross section may be smaller, because 
the incident momentum is considered to be small and 
the process will be near threshold. From this estimate, 
if we take into account the broadening of the width, the 
mass shift becomes slightly smaller, by about 0.5 MeV, 
according to the results shown in Fig. [151 Therefore, this 
justifies the argument in Ref. [2l| that the influence of 
the decay widths is expected to be small. 

The change of spectral properties in the nuclear mat- 
ter can be experimentally investigated by Panda experi- 
ment at GSI-FAIR in which incident anti-proton collide 
with nuclear target. Here we present some predictions 
for cross sections of charmonium production through pp 
annihilation and subsequent decay into dileptons or ra- 
diative decay in the experiment. We compute the cross 
sections with the Breit-Wigner formula 



cbw(s) = 



B in B out (2J+l) 4tt 



sT 2 



(2 Sl + l)(2s 2 + 1) fc c 2 m (s - m 2 ) 2 + S r 2 ot+mcd : 



(31) 

where s, fc 2 m and m are the Mandelstam variable, cm. 
momentum and mass of charmonium with spin J, re- 
spectively. r to t is the total decay width of the charmo- 
nium and Ftot+med = Ttot "i" F mec [ium- 

.Bin and B out are 

the branching fraction of the resonance into the entrance 
and exit channels. Sj is the spin of the incident parti- 
cles, which are anti-protons and protons in the present 
calculation. Since the target protons are in nucleus, we 
have to take the Fermi motion into account for accurate 
estimation. We average the Breit-Wigner cross section 
with respect to target momentum as 



Po Jq 



- k 2 



, dkdtt 



(27T) 



rCBW- 



(32) 



In addition to J/ip and r) c , we also calculate cross sec- 
tions for Xc which are expected to show larger mass shift 
8m ~ —40 ~ — 60 MeV [51j. Parameters in the calcula- 
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TABLE III: Parameters and results in charmonium productions at GSI-FAIR. Cross sections and event per day correspond to 



the case of 


maximum medium 


width, r mcd = 


20 MeV. 








riesonance 


771 [ivie V J 


o//2[ivie v j 


1 tot 


r lllcLl OhQXG 


o"bw stt peak 


Events per day 


J/V 


3097 


-7 


93.4keV 


e + + e~~ 


0.435 pb 


7.5 


Vc 


2980 


-4 


25.5MeV 


e+ + e~ 


10.7 pb 


184 


XcO 


3415 


-60 


10.4MeV 


J/ip + j 


18.0 pb 


311 


Xcl 


3511 


-60 


0.89MeV 


J/ip + j 


4.5 pb 


78 


Xc2 


3556 


-60 


2.05MeV 




19.8 pb 


343 



pp -> JA|/ 



eV 



(5m,r med )=(0,1) 
[MeV] =(0,5) 
=(-7,1) 
=(-7,5) 




(5m,r med )=(0,10) 
[MeV] =(0,20) 
=(-7,10) 
=(-7,20) 
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FIG. 16: (Color online) Cross section of J/ip production in 
p — A collisions. Upper panel shows smaller medium width 
(1 and 5 MeV) cases and lower one shows larger (10 and 20 
MeV) cases for with and without mass shift. 



tions are summarized in Table HTT1 r mc di U m is treated as 
a free parameter varied from 1 MeV to 20 MeV. 

Results of the cross sections as a function of incident 
anti- proton energy are shown in Figs. I16H20I We can 
clearly see that sharp peaks of the resonances disappear. 
This is because of the Fermi motion of the target pro- 
tons in the nucleus. For example, incident energy to 
create J ftp (3097) is E lab = 4.17 GeV, but the fluctua- 
tion of the target momentum makes it possible to create 
J/ip with 3.17 < E lah < 5.51 GeV, in which the min- 
imum and the maximum E^b correspond to the target 
momentum along the collision axis p2 Z = —k-p and fcp, 



FIG. 17: (Color online) Same as Fig. US] but for 7] c . 



respectively. This effect considerably broadens the cross 
section. Consequently, one needs no fine tuning of in- 
cident proton energy to produce charmonium. For J/ijj 
and rye, mass shifts are so small that the peak positions 
of incident energy do not change. However, mass shift of 
Xc, ~ —60 MeV, is sufficiently large to show clear shift 
of the peak in the cross section. In these calculations, we 
treat r mc( jium as a free parameter. It is shown that this 
parameter affects only on the magnitude of the cross sec- 
tion, which is larger for smaller change from the vacuum 
width. Hence, though we cannot predict both of mass 
shift and in-medium width, we can obtain information on 
both quantities from the experimentally measured cross 
sections. We summarized the cross sections and expected 
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FIG. 18: (Color online) Same as Fig. [TS] but for XcO- 



FIG. 19: (Color online) Same as Fig. [16] but for Xd- 



event rate at GSI-FAIR, of which luminosity is expected 
to be 2 x 10 32 cm _2 s _1 , in last two column of Tabic. HTT1 
We can see that the expected event rates are large enough 
for the mass shift of Xc to be observed. 

Finally we address possible improvements of this work. 
Since we restricted ourselves to the hot medium which 
consists of gluons only in the first part of this paper, we 
should take the quarks into account for more realistic 
estimation. 

To consider the quark effects, first we consider the 
quark operators appearing in the OPE side. We can ne- 
glect the light quark contribution to the OPE, because 
the light quark operators appear in the OPE at order 
a 2 s (q 2 ): This is why the light quark condensate can be 
neglected in the sum rules for heavy quark system. On 
the other hand, thermal heavy quarks that directly cou- 
ple to the heavy quark current contribute to the OPE 
at leading order. This is different from the heavy quark 
condensates that are perturbatively generated from the 
gluon condensates, and contribute to the OPE through 
gluon condensates, whose Wilson coefficients are calcu- 
lated in the momentum representation. The direct ther- 
mal quark contributions are called the scattering terms. 
However, similar terms also appear in the phenomenolog- 
ical side, which also has free charm quark mode that is 
not coupled with a light quark in the form of a D meson 



above T c as been recently studied in Ref. [52j. Therefore, 
the scattering term will cancel out between the OPE side 
and the phenomenological side in the deconfined medium. 

Second, the gluon condensates themselves can have a 
different temperature dependence in the presence of dy- 
namical quarks. As discussed before, the important in- 
put for the mass and width change is the temperature 
dependence of gluon condensates in Fig. [2j in particular 
the dominant contribution comes from the temperature 
dependence of Go- For that purpose, we note that the 
trace of the energy momentum tensor to leading order is 
given as, 



11 - 2/3ATA 
8 J 



^m q (qq). 

(33) 

Therefore, we start from the lattice calculation of the 
trace of the energy momentum tensor for the full QCD 
with realistic quark masses given in Ref. [53j. Then, we 
subtract the fermionic part of the trace anomaly, which 
was also shown in the literature, from the total. Next, we 
divide the result for the relevant prefactor with Nf = 3 
multiplying the gluon condensate as given in Eq. (|33|) . 
Since the critical temperature T c differs, we compared it 
as a function of T/T c in which T c = 196 MeV for the 
full QCD case [53| . As can be seen in Fig. [5TJ the magni- 



16 



JA|/y 



n 

Q. 

t> 

c 
o 

o 

CD 

CO 
CO 

o 



40 
30 
20 
10 




(5m,r med =(0,1) 
[MeV] =(0,5) 
=(-60,1) 
=(-60,5) 




(5m,r me J)=(0,10^ 
[MeV] =(0,20) 
=(-60,10) 
=(-60,20) 




■lab 



[GeV] 



FIG. 20: (Color online) Same as Fig. [16] but for Xc2. 



tude of the resulting change near the critical temperature 
are remarkably similar between the full and pure gluon 
QCD; although the slope at T c is milder for full QCD 
as a consequence of rapid cross over transition instead 
of the first order phase transition. Since the change of 
the condensate sets in at a lower T/T c in the full QCD 
case, the mass and width of charmonia might start vary- 
ing at a lower temperature in the realistic case than in 
the pure glue theory. As for the twist-2 condensates, re- 
sults will not be affected so much by taking into account 
the fermionic part since the effect will be small at this 
temperature region. Therefore we believe our main ar- 
gument and the quantitative result will not be alter even 
in the realistic situation. 

It is also important to study change of Xc at finite tem- 
perature, which may influence the quantitative feature 
of the sequential melting (54[, since non- negligible frac- 
tion of J/ip comes from decay of ip' and \c hi relativistic 
heavy ion collisions. This can be done by calculating 
Wilson coefficients for tensor operators for these chan- 
nels. It should be also noted that the continuum part of 
the spectral function may play an important role. This 
will be possible by modeling the medium with a gas of 
quasi-particle. One more thing to be done is the exten- 
sion to higher temperature. The failure of T > 1.06T C for 
J/ip and T > 1.04T C for r\ c originates from the instability 
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FIG. 21: (Color online) Comparison of the scalar gluon con- 
densate in the pure gauge theory with the one of full QCD. 
Horizontal errorbars in the full QCD case are drawn by as- 
suming the 2% uncertainty in the conversion from the lattice 
units to physical temperature [53 ]. 



of the moment ratio of the OPE side. The twist-2 gluon 
condensates becomes larger as temperature increases and 
then leads to the breakdown of the stability in the OPE 
side including up to dimension 4 (see Fig. O and 0(a s ). 
In Fig. [Jj we can also see that the expansion is not good 
at large n that stabilize the moment ratio at higher tem- 
perature. These facts suggest the necessity of includ- 
ing higher dimensional operators, which is examined in 
Ref. [551 ] - However, we do not know a simple way to 
extract the temperature dependencies of the higher di- 
mensional operators from the lattice calculation, as was 
done in the present work for dimension 4 operators. The 
other way of the extension is to improve the phenomeno- 
logical side such that it includes temperature dependent 
continuum contribution. The decrease of the scalar gluon 
condensates above T c indicates perturbative contribution 
becomes more important at higher temperature. If we 
can construct a more appropriate phenomenological side 
reflecting the nature of the strongly interacting matter, it 
will lead to n-independent results for physical parameters 
until the charmonia really dissolve. 

We also note that there are some spaces to improve 
the analyses for nuclear matter. Especially, the present 
analysis shows the mass shift of \ c states are likely to be 
observed in the forthcoming experiment. However, the 
current estimate of the mass shift is not a decisive one; 
we have to take the twist-2 contribution into account for 
a more accurate estimation. 

In summary, we have given a comprehensive analysis 
on medium-induced change of the spectral properties of 
J /if) and i] c in the hot gluonic medium and the nuclear 
medium by making use of QCD sum rules. In the case 
of the gluonic medium, our analysis shows there must 
be a notable change of mass or width, or both around 
T c , caused by the rapid change of the gluon condensates. 
Although the present formalism is found to be applicable 
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only up to T ~ 1.06T C , the change of mass and width can 
maximally reach a few hundred MeV. We have discussed 
its implication for future heavy ion experiment at CERN- 
LHC. As for the nuclear matter case, we extend the past 
works to include small but finite width and check the 
robustness by varying the scale parameter of the theory. 
We also examined the possibility of detecting such mass 
shifts in the future experiment at GSI-FAIR. Although 
J/ip and r\ c do not show prominent signals, Xc exhibits 
more promising results. These analyses give the basis of 
future improvements to study the nature of the strongly 
interacting matter deeply with charmonia. 



the Korean Research Foundation KRF-2006-C00011 and 
by the Yonsei University research fund. K. M. would like 
to thank the members of the high energy physics group 
of Waseda University for allowing him to use their work- 
stations. We also would like to acknowledge T. Hatsuda 
for his fruitful comments and discussions. 



APPENDIX A: WILSON COEFFICIENTS 



Acknowledgments 

This work was supported by BK21 Program of the Ko- 
rean Ministry of Education. S. H. L. was supported by 



Here we list explicit forms of the Wilson coefficients 
which appear in Eq. (|10|) and are originally given in 
Refs. [l| and [Hj]. In the following, p = £/(l + £) and 
F(a, b, c; x) is the hypergeometric function 2 Fi(a, b, c; x). 



J 



For the pseudoscalar channel, 
3 2"(n-l)! 



<(0 



8tt 2 (2n+l)!! 
(2n + l)!! 



(4m 2 )-"(l + tr n F(n, 1/2, n + 3/2; p) 
1 



3 • 2 n " 1 n! L" 2(n+ 1) V2 4,t 
1 |"8 4 F(n,3/2,n + 3/2;p) 
3 



-tt--^- )F(n,l,n + 2;p) 



1 



l/2,n + 3/2;/?) 
1 F(n, 3/2, n+ 5/2; p) 



1 3 

— I — 7T 

2 4tt 



2n 



n F{n, 1/2, n + 3/2; p) 6n + 3/2 F(n, 1/2, n + 3/2; p) 
ln(2 + Q (2 + Q F(n + 1, 1/2, rc + 3/2; p) 
7T (1+0 2 F(n,l/2,n+3/2;p) ' 



n(n + l)(n + 2)(n + 3) 



2n + 3 



F(n + 1, -3/2, n+ 5/2; p) 6 F(n + 1, -1/2, n + 5/2; p) 



F(n, 1/2, n + 3/2; p) 



F(n,l/2,n + 3/2;p) 



<£(0 =&rft) 



4n(n + 1) F(n + 1,-1/2, n + 3/2; p) 



(1+0 F(n, 1/2, n + 3/2; p) 
Similarly, for the vector channel, 

3 2 n (n + l)(n-l)!F(n, l/2,n + 5/2; p) 



&Jft) 



"4tt 2 



(2n + 3)!! 
(2n + l)!! 



[(4m 2 )(l + 0] n 
'2n + 3 



3 • 2 n ~ 1 nlF(n, l/2,n + 5/2; p) V2n + 2 



+ 



F(n,2,n + 3;p) f n 3 
3(n+ l)(n + 2) V2 ~ 4^ 











7T 










3 \ 






- 


- 2n 


47T I 





1 /7T 



2 V 2 



3 

47T 



F(re,l,n + 2;p) 

71 + 1 



ln(2 + Q (2 + Q Fjn + 1, 1/2, n+ 7/2; p) 
7T (1 + 2 F(77,l/2,77 + 5/2;p) : 



i(n + !)(»+ 2)(n + 3) F(n + 2, -1/2, ti + 7/2; p) 



(2n + 5)(l + 2 F(n, 1/2, n + 5/2; p) 

v (f] =h v (n _ 4n(n+l) Fjn + 2, 3/2, n + 7/2; p) 
c„UJ „W 3(2n + 5)(l + 2 F(n, 1/2, 7i + 5/2; p) " 

In Eqs. (|A1|) and (|A5|) . m is the running quark mass m = m c {p 2 



-(£ + l)7n 2 ) which is given by [56 



™ C (Q _ 1 _ Og 
TO c(£ = 0) 7T 



2 + e 



ln(2 + 0-21n2 



(Al) 



(A2) 
(A3) 
(A4) 
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(A9) 
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